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We propose a general framework to solve tight binding models in D dimensional lattices driven 
by ac electric fields. Our method is valid for arbitrary driving regimes and allows to obtain effective 
Hamiltonians for different external fields configurations. We establish an equivalence with time 
independent lattices in D+l dimensions, and analyze their topological properties. Further, we 
demonstrate that non-adiabaticity drives a transition from topological invariants defined in D+l 
to D dimensions. Our approach provides a theoretical framework to analyze ac driven systems, 
with potential applications in topological states of matter, and non-adiabatic topological quantum 
computation, predicting novel outcomes for future experiments. 
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Introduction: Periodically driven quantum systems 
has been a fastly growing research field in the last years. 
The development of effective Hamiltonians describing ac 
driven systems at certain regimes, has allowed to pre- 
dict novel properties such as topological phases [IH3, and 
quantum phase transitions^, 5] that otherwise, would be 
impossible to achieve in the undriven case. Therefore, 
the application of ac fields has become a very promising 
tool to engineer quantum systems. On the other hand, 
the obtention of effective Hamiltonians can be a difficult 
task, depending on the driving regime to be considered. 

In this work, we provide a general framework to study 
periodically driven quantum lattices. By means of our 
approach, it is possible to solve with arbitrary accuracy 
their time evolution, and obtain effective Hamiltonians 
for the different driving regimes. We consider solutions 
of Floquet-Bloch form, based on the symmetries of the 
system, and characterize the states in terms of the quasi- 
energies, which are well defined for all driving regimes. 
The states belong to a composed Hilbert space, in which 
time is treated as a parameter. As we will see below, 
it allows to formally describe the ac driven D dimen- 
sional lattice, as analog to a time independent D+l di- 
mensional lattice. This description enlightens the under- 
laying structure of periodically driven systems, in which 
the initial Bloch band splits into several copies (Floquet- 
Bloch bands), where the coupling between them directly 
depends on the driving regime. Interestingly, the isolated 
Floquet-Bloch bands possess the same topological prop- 
erties as the Bloch bands of the undriven system, now 
tuned by the external field parameters. Thus, the topo- 
logical invariants for the isolated Floquet-Bloch bands 
can be obtained following the general classification of 
time independent systems (AZ classes However, 
we also demonstrate that lowering the frequency, the 
bands couple to each other. In that case, the topological 
structures are classified according to a base manifold of 
dimension D+l. In consequence, one can simulate higher 
dimensional tight binding (TB) models with exotic tun- 
able hoppings by just coupling the system to ac electric 
fields. 

Our approach is valid for arbitrary dimension, and it 



allows to independently analyze the effect of the field am- 
plitude and frequency. We show that the field amplitude 
controls the renormalization of the system parameters, 
while the frequency acts analogously to a DC electric field 
in the extra dimension. In particular, this last property 
relates the high frequency regime with the existence of 
Bloch oscillations and Landau-Zener transitions between 
bands|^|, establishing a direct relation between diabatic 
regime and localization[10(. We illustrate our formalism 
with the analysis of an ac driven dimers chain 11 - 13| ■ 
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Figure 1: Schematic figure of the equivalence between a peri- 
odically driven ID chain, and the effective time independent 
lattice in 2D. We plot the positions of the undressed states 
(blue) and their spatial distribution at sites ra, r G Z. Each 
site is coupled to a set of dressed states (red color), with 
coupling oc tJ p (Ao), being p = n — m the difference in the 
number of photons between the coupled Floquet bands, and 
Ao the vector potential amplitude. We draw the Wigner-Seitz 
unit cell in 2D (green), and the effective dc electric field (see 
text) along the energy axis £ (green dotted arrow). The extra 
dimension £ arises due to the time periodicity. 

Theory: We start by considering a Hamilto- 
nian with lattice and time translation invariance: 
H (x + a. t ,t + T) = H (x + ai,t) = H (x, t + T), charac- 
terized by lattice vectors a.; and time period T = 2ir /uj. 
Under these assumptions, we assume the Floquet-Bloch 
ansatz0, El |* a ,k(x,t)) = e lk x —° k > Q , k (x, t)), 
being e Qj k the quasi-energy for the a Floquet state, and 
k the wave-vector. The Floquet-Bloch states \u a ,k (x, t)) 
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are periodic in both x and t, and belong to the composed 
Hilbert space (Sambe space) S = Jf? (8) T[lB|, where 
T is the space of T-periodic functions and the 
Hilbert space. It defines a composed scalar product 
((. . .)) = Jq T (- • .)dt/T. In this space, the Floquet states 
fulfill the orthogonality condition ((«„ (t) \up (t))) = 
T^n.m Jo elu(n ' m)t ( u a,n\ u /3,m)dt = 8 a ^5 n , m , and the 
Floquet equation: 

H (k,t) |u a ,k) = ea,k|Ua,k), (1) 

W(k,t) = e" ik ' x (if (t) - idt) e* x 

= H k (t)-id t (2) 

where |w a ,n) are the coefficients of the Fourier expansion 
in t, and W(k,t) the Floquet operator. We also define 
the Floquet-Bloch annihilation and creation operators 
| c k ' c Wl k}- Note that the Floquet ansatz allows 
to map the time dependent Schrodinger equation to an 
eigenvalue equation, in which t and k are both parame- 
ters. We propose below that an exact mapping between 
a D dimensional ac driven system and a D+l undriven 
one can be established, in which H (k, f) is equivalent 
to an static Hamiltonian. It will allow to classify the 
topological invariants in terms of the mappings from the 
parameter space to the set of Floquet operators H (k, t) , 
being the parameter space now given by S 1 x P, where 
n is the dimension of the First Brillouin zone (FBZ). 
We first define the Fourier transforms: 

N oo 

c(*) a ,k = N- D / 2 J2 E ^>e <k ' R ' +W (3) 

j—1 n— — oo 

N oo 

c(t)U = n- d ^ E 

j=l n— — oc 

where N is the number of sites in the the lattice with 
periodic boundary conditions and D the dimension of 
the undriven system. For a time dependent Hamilto- 
nian in the dipolar approximation, the relation between 
undriven and driven system is given by the minimal cou- 
pling k — > K (t) = k + A (t) , where A (t) is the vector 
potential. Equivalently, one can notice that in the un- 
driven TB Hamiltonian, the hopping parameters Tjj are 
k independent. Thus, one can obtain the time depen- 
dent TB by using the minimal coupling in the inverse 
Fourier transform, leading to the time dependent hop- 
pings Tjj (t) = Tjje iA ^'^ R i^ Rl ' , and being the Hamilto- 
nian 

ff k(t)=EE T ^(*) e<k ' (R3 '" Ri)c (*)«,k C (*)a ) k, (4) 
a,k j,Z 

where the a index labels the Floquet state, and the time 
dependence in the operators is included. Including the 



expansion in Fourier series of c (t) a k , and c (t) k (EqJ3|) 

^(*) = EEE T ^ (*) ik.nw^^^-^ 

a,k n.m jj 

(5) 

being the quadrivectors k (t) = {—t, k) and p n j = 
(nu), Rj). EqJS] gives a description of the time dependent 
Hamiltonian in terms of the time independent operators 
|c ai k,n,c^ kn |. Finally, the use of the composed scalar 
product allows to obtain the quasi-energies by diagonal- 
ization of the matrix: 

((u a ,k,n\H (k, t) |u Q ,k,m)) = f n ,m ~ nu5 n , m , (6) 

i r T N 

f n , m = - / Y, (t) e^-^-^dt, (7) 

where nuj8 n ^ m is the Fourier space representation of — idt- 
EqlHl is analog to a time independent TB in D+l di- 
mensions with an electric field of unit intensity applied 
along the extra dimension, and sites labeled by (n, j) (see 
FigHJ. The effective electric field breaks translational 
symmetry in the £ axis, and one can differentiate two 
regimes: Low and high frequency. 

In the low frequency regime (w <C Tjj), it is a good 
approximation to neglect the effect of the time derivative 
in the Floquet operator (Eq[2]) , or equivalently, the Stark 
shift due to the effective electric field. Then, we restore 
the £ axis translational symmetry, and the operator can 
be diagonalized by Fourier transform of Eq(B]to t domain. 
The obtained T-L (k, t) is analog to a Hamiltonian over a 
D+l compact base manifold, that we define as the First 
Floquet Brillouin Zone (FFBZ), parametrized by {t, k} £ 
S 1 x T™ [l3|. Hence, the system topology is classified 
according to the AZ class of D+l static Hamiltonians. 
Note that if one initially assumes adiabatic evolution, 
the Floquet structure, the parameters renormalization by 
the field amplitude, and the additional dimension of the 
base manifold are not obtained. Thus, the topological 
classification is not clearly stablished. 

For high frequency (cj S> Tjj), the effective elec- 
tric field produces Bloch-Zener transitions and Bloch 
oscillations 0, [Io| . inducing localization in raj, and de- 
coupling the Floquet bands. Then, Eq|6] becomes block- 
diagonal in Fourier space, and the effective Floquet op- 
erator time independent. In addition, it is defined over 
a D dimensional base manifold k E T™ (FBZ). In conse- 
quence, the topological classification is given by the AZ 
classes for time independent systems in D dimensions. 
The transition from a D+l to a D dimensional base man- 
ifold (equivalently, from the FFBZ to the FBZ), as one 
increases the frequency uj, is driven by the non-adiabatic 
processes along the t axis[l8|. 

Importantly, note that in general the hopping between 
n, m neighbors depends on the amplitude of the vector 
potential (Eq[7]). It allows to design "effective lattices" 
by tuning the hoppings with the ac ficld[19]. 
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Periodically driven dimers chain: Here we consider 
the case of a dimers chain coupled to an ac electric 
field, with hoppings t and r', and periodic boundary 
conditions (FigHJ.The electric field E (t) = -d t A(t), is 
given by the vector potential A(t) = Ao sin (ujt), where 
Aq = qEo/uj (we fix q — — 1). 
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Figure 2: Schematic figure for a dimers chain with two in- 
equivalent atoms (A, B) for unit cell (green area), bo is the 
distance between A-B atoms within the same primitive cell, 
ao the lattice translation vector, r' the hopping within the 
same dimer, and r the hopping to the next one. 



Figure 3: Quasi-energy spectrum vs Eq/uj for A = 0.3, and 
bo — 0, considering Eq(9] (high frequency). The band struc- 
ture is obtained by considering 10 k vectors equally spaced 
within the FBZ. We also included in green color the gapless 
modes obtained from the numerical calculation, ui — 10. All 
parameters in units of r. 



By calculating f n>m from EqJTl one obtains: 

r. =r( ° PF{k) ) (8) 
Tn ' m ~ \~p F (k) J ' W 

p F (k) = \e- lkb «J n _ m (A b a ) 

+e ik(a -b )j m _ n {Aq ( ao _6 o))) 

p F (k) = Ae ifcb °J m _„(A & ) 

+e - ik{ao - bo)jn _ m {Aq ( ao _ 6o)) . 

where A = r'/r, and (n,m) £ Z. In contrast with the 
undriven case, the spectrum depends on the intra-dimer 
distance bo, and the hoppings are renormalized by the 
field amplitude. Note that for the limit of vanishing ac- 
field Ao — > 0, the quasi-energies match the energies of 
the undriven system. 

We first consider the high frequency regime (oi 
t, t'), and select the Floquet band n = m = 0. In this 
case, the Floquet operator is block diagonal, and the sys- 
tem can be described by a time independent 2 by 2 matrix 

n[ 0) =rg(k)-a, (9) 

where g(k) = (Si (pf) , 3 (Pf) , 0) for n = m = 0, and 
& = {°~x, o~ y , a z ) are the Pauli matrices. 

In order to compare with the case with periodic bound- 
ary conditions discussed above, we also consider a finite 
dimers chain with Hamiltonian H (t) — Ho + qE (t) x, 
being Hq the time independent TB Hamiltonian, and 
qE (t) x the coupling with the electric field (details in 

m. 

FigOS shows the quasi-energy spectrum in high fre- 
quency regime. We included in green dotted lines the 
regions of existence of edge states, obtained from the nu- 
merical calculation of the finite size system. 

The appearance of zero energy modes is a finite size 
effect linked to the underlaying topology of the sys- 
tem. The relation between them and the bulk topol- 
ogy is via the bulk to edge correspondence, which relates 



the number of zero energy modes at the boundary, car- 
rying a topological number, with the bulk topological 
invariants |20j. Therefore, the calculation of the topo- 
logical invariants of %^ should predict their existence 
in this regime (FigJ3]). Wi belongs to the BDI class 
(as the one corresponding to the undriven system), with 
time reversal, particle-hole, and chiral symmetry [8j. In 
ID, the BDI class is characterized by a winding number 
v\, which classifies mappings T 1 — > R 2 — {0} ~ S 1 , from 
the FBZ to the family of Hamiltonians : 

vi = <£{u at k\idk\u a ,k)dk 

= |(l + S ign(J 2 (t/)-A 2 J 2 ^))), (10) 

where y = Aq (ao — bo), x = Aobo, and |u a .fe) are the 
closed lifts of H^] ■ Eq(Tn]shows, that in contrast with the 
undriven case[ll). one can create non-trivial topological 
phases even for A > 1, where the undriven system is 
in the trivial phase (FigQ] left) . This is an example of 
topology induced by the driving. 

In Fig0] (right) we plot the phase diagram correspond- 
ing to A = 0.3, which correctly predicts the existence of 
edge states for bo — (Fig(3|). 

In summary, we have shown that in the high frequency 
regime, the topological properties can be obtained using 
an effective static Hamiltonian ni , and that they can 
be tuned by the field amplitude. 

As we decrease ui, the different Floquet bands couple to 
each other, and the isolated band picture is not accurate. 
In this regime, one must consider the full Floquet opera- 
tor (Eq[T]), which for this system is not exactly solvable. 
Due to the coupling between Floquet bands, two different 
but related effects happen as lo is reduced: Bands inver- 
sions, and the emergence of a D+l parameter space. 

Bands inversions correspond to crossings of the bands, 
in which the symmetry is exchanged (e.g., it occurs 
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Figure 4: Topological phase diagram in the high frequency 
regime for an ac driven dimers chain. We considered A = 1.5 
(left), and A = 0.3 (right). Dark color area means v\ — ir, 
and light area means v\ = 0. Note that even for A > 1 we can 
induce a non-trivial topology, in contrast with the undriven 
case. Further, the phase diagram for A = 0.3 (right) agrees 
with the existence of edge states in FigO (bo = 0). 



Note that in difference with EqJHl g (k, t) now depends 
on t, and v\ can be defined along the two inequivalent 
axis of the torus 

1 f d 



v\ (v) 



— 7T 
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in quantum wells of HgTe/CdTe as the well thickness 
reaches a critical value|20]). This effect happens in ac 
driven systems as we decrease the frequency, because the 
distance between Floquet bands is proportional to w. If 
the maximum width of an isolated Floquet band, is given 
by Se < u) (e a £ [— uj/2, w/2]). Then, the Floquet bands 
at ±u> close the gap when cj = Se/2. As a general rule, 
band inversions happen for every: 



UJ. n 



ii 

2n' 



n e 1 + , 



(11) 



where Z + denotes the set of positive integers. Therefore, 
at u n the ±nu Floquet bands close the gap, switching 
between an ordinary and a topological insulating phase. 
For example, for a dimers chain with bo = the max- 
imum width coincides with the undriven system band 
width, i.e., Se = SE = 2t\/1 + A 2 + 2A. Then, by means 
of Eq llll it is possible to track the bands inversions in 
terms of the undriven system[l7j • 

For uj -C r, r', a large number of bands inversions oc- 
cur, being difficult to track all of them. In addition, the 
presence of a D+l base manifold becomes important. In 
that case, one can neglect the time derivative in Eqf2J and 
diagonalize the operator in t domain. In that case, one 
obtains a 2 by 2 Floquet operator H (k,t) ~ H (k,t) NN 
defined over the FFBZ, where H (k,t) NN extends up to 
next nearest neighbors coupling in (n, m) (it is a good 
approximation for Aq < 1). However, 2D Hamiltonians 
in the BDI class are topologically trivial, and only the ID 
topological invariant V\ is still non zero. It means that 
all changes in the topological properties will be reflected 
in v\. In addition, for BDI Hamiltonians one can com- 
pute the winding number gra phicallyjllj, in terms of the 
divergences of the phase <fi(k,t) = arctan (g y /g x ) in the 
FFBZ (Figl5]), being g XtV the components of the vector 



H(k,t] 



NN 



rg(k,t) 



(12) 



Figure 5: Plot of (k, t) all over the FFBZ for A = (left) 
and Aq — 2 (right). Paths parallel to k cross two discontinu- 
ities meaning V\ = ±7r (red arrow), on the other hand, paths 
parallel to t have v\ — because they wind back and forth 
(black arrow). 

For the case of a finite chain, the existence of bound- 
ary states depends on the loops along the k axis, i.e., on 
v\ (i)[U. In Figj5]^i (t) = tt for all t, independently of 
the value of Eq/oj (trajectories always cross two discon- 
tinuities). Thus, gap inversions are not relevant for the 
existence of edge states in the low frequency regime, be- 
ing a feature of the transition from the FBZ to the FFBZ. 
This can be seen numerically in the finite size system, in 
which zero energy modes are present independently on 
the gap inversions [l7) . 

Conclusions: We have derived a general approach to 
solve periodically driven D dimensional lattices. It allows 
to obtain effective Hamiltonians for the different driving 
regimes and a complete topological classification in terms 
of AZ classes. We show that even for high frequency, the 
underlying topology of the undriven system is present due 
to the time periodicity. In addition, we show that at low 
frequency, the Floquet operator is analog to the one of a 
static system in D+l dimensions, leading to interesting 
topological states of matter which otherwise would be 
inaccessible. Finally, we also described the mechanism 
of bands inversion in ac driven systems and its relation 
with the topology of the system. 

We have shown a large horizon of possibilities for peri- 
odically driven systems, which in addition, can be studied 
experimentally e.g., by measuring the electric polarizabil- 
ity or the appearance of boundary states. The driving al- 
lows to simulate properties of undriven systems in higher 
dimensions and the obtention of new topological phases 
due to tunable hoppings[19j. One could also think in 
more exotic types of zero energy modes in the low fre- 
quency regime, as for example those in the boundary 
between a driven and undriven materials. Also, further 
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studies in the case of bichromatic ac fields could be inter- 
esting. Finally, the results in the high frequency regime 
can be used for non-adiabatic topological quantum com- 
putation. 



I. SUPLEMENTARY INFORMATION 

S.l: Relation between undriven and ac driven 
Hamiltonians 
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For a periodic lattice, driven by an ac electric field 
under the dipolar approximation (A (x, t) ~ A(t)), the 
relation between the static Hamiltonian and the ac driven 
Hamiltonian can be obtained by noticing that: 

,2 



H 
H k 

H(t) 
flk(t) 



-ik.-x.Tjik.-x 



(13) 



(p+Awr 



(P±k) 

2m 



+ V(x), 



2m 



+ V(x), 



-ik-x 



H(t)e 



ik-x 



(p + k + A(<)) 2 



2 m 



_ e -i(k+A(t))-xjj e t(k+A(i)).x_ 



It relies in the minimal coupling k — > K (t) = k + A (t). 
Explicitly for the case of a tight binding Hamiltonian, 
one should notice that the hoppings can be defined as 
the Fourier transform of the energy: 



= iV- 1 ^ £ k e- 4k - (R ^ Ri) 
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which are k independent. Therefore, the inverse Fourier 
transform of Tjj can encode the time dependence in the 
minimal coupling as: 



H k (t) 



^ D E^ elK(tHR ^ R,)c wL k c W«, k ' 

3,1 



where the time dependent hoppings Tj t i (t) = 
Tj t ie lA ^'( R3 ~ R ') are obtained from the time inde- 
pendent ones. This result relies on the same principle as 
the Bloch equation for a particle in a dc electric field [1]. 

S.2: ID system in low frequency regime 

Let us consider as an example, a ID chain driven by an 
ac electric field, where u> <C tjj, being T\j the hopping 
between sites, and uj the frequency of the ac field. In 
this case, one can neglect the time derivative term of 
the Floquet operator, or equivalently the effective static 
electric field, as we described in the main text. Then, 
the translational symmetry along the energy axis E is 
recovered, and one can write the Floquet operator in t 
domain, where it is diagonal: 

a,k m,n 

being M a normalization factor. Thus, because the Flo- 
quet operator (Eq |14p depends on the compact parame- 
ter space (fc, t), we can classify the topological properties 
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of % (k, t) according to the AZ classification of time in- 
dependent D+l dimensional systems [2j. It classifies the 
mappings from the torus to the family of Floquet oper- 
ators W (k,t). For example, the case of a 2 band model 
with Floquet operator: 

H (k,t) = h{k,t) ■ a, 

possess a first Chern number given by: 



''' = / / h(k,t) ■ ( -^-h(k,t) x J^MM) ) dtdk 



dk 



dt 



d - d ~ 

h (k, t) ■ ( — h (k, t) x — h (k, t) ) dtdk, 

FFBZ V &k Ot 



being the First Floquet Brillouin Zone (FFBZ) homeo- 
morphic to a torus T 1+1 , and h = h/\h\. Note that for a 
one dimensional time independent system, c\ would al- 
ways vanish because the FBZ is given by T . However, 
the increase of the parameter space dimension, due to 
the time periodicity, allows for higher order topological 
invariants. 

In the particular case of a dimers chain, discussed in 
the main text, c\ — 0. The reason is that the Floquet 
operator belongs to the BDI class in 2D. Therefore, the 
winding number v\ is the topological invariant that dif- 
ferentiates our system from an ordinary insulator. 



S.3 Dimers chain 
S.3.1: Hamiltonians and topological invariants 

The undriven tight binding model for nearest neighbors 
is given by: 



p(k) 
p(k)* 

p{k) = T 'e- ikho + T e lk{a °- bo \ 



(15) 



E± = ±ry/\ 2 + 1 + 2Acos (ka ), 

being A = r'/r (see Fig. 2 in the main text). Importantly, 
the dispersion relation does not depend on 6q, and for the 



condition bo = ao/2, and r = r' in one recovers the 
energy spectrum of the linear chain. 

We consider the ac vector potential A (t) = Aq sin (u>t), 
being Aq = qEo/uj. By means of the minimal coupling 
we arrive at the time dependent Hamiltonian: 



H K(t) 



p(k,t) 



w - ■ y P (k,ty o 

p(k t) = Xe~ *(*+ J *° sin(ojt))6 _|_ e i(k+A sin(w*))(oa-6o) 

In order to calculate f„, m (Eq.7 in the main text), we use 
the identity J n (x) = ± ^ e ix*H"t)-iput d ^ i ea ding to: 

Tn > m = T [p F (k) J' (16) 
p F {k) = \e- lkbo J n - m (A b ) (17) 
+e iHao-b )j m _ n ( Ao ( ao _6 )), 

p F (k) ee A e ' ifcbo J m _„(A 6o) 

+e - ik ( ao - bo)jn _ m (Aa ( ao _6 o)) . 

Finally, one can obtain the matrix elements of the Flo- 
quet operator as: 



o/(n,m) _ - 

n k — t„ 



nuj5 n , m , 



(18) 



where we have included the time derivative operator in 
Fourier space. Note that this is an infinite matrix because 
(n, m) € Z. In the high frequency regime (w 3> r, r'), the 
second term in the right hand side of EqJTH] dominates, 
and the matrix is approximately block diagonal. Thus, 
we select the Floquet band m — n = for simplicity, be- 
ing the effective Hamiltonian given by the 2 by 2 matrix: 



(0) 




H 



pf ee \J (A b ) + e lka °J (Ao(ao-bo)). (19) 

Note that we have considered p^p in a different ba- 
sis than in EqJTTl Both basis differ in a phase factor 
e ±ikb _ rj^g reagon j s i\y a i i n order to properly obtain 

the topological properties we must consider the closed 
lifts basis[3|, |J|. The quasi-energy spectrum at high fre- 
quency is (q = —e = — 1): 



e ±,fc 



±T 



J o (y) + 2 ^ cos ( ka o) Jo (x) J (y), 



(20) 



where x ee ^o^Oi an d y ee ^4o(ao~^o)- Fig [S] shows 
a comparison between the numerical calculation of a fi- 
nite size system in high frequency regime (right), and 



the quasi-energies obtained in the high frequency regime 
for the model with periodic boundary conditions, Eq |20l 
(left). 
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Figure 6: Quasi-energy spectrum vs Eq/oj for A = 0.3, and bo = 0, considering Eq |20l for the model with periodic boundary 
conditions (left), and the numerical solution for the finite system with 20 dimers (right). For the finite dimers chain we have 
considered n,m — 150 sidebands and uu — 10. Note the existence of gapless modes for a certain range of Eq/uj. All parameters 
in units of r = 1. 



Out of the high frequency regime the Floquet bands cou- 
ple, and we must solve the full Floquet equation. Fig|7] 
shows a comparison between the model with periodic 
boundary conditions (left) and the finite one (right), for 
frequency u ~ r, r'. In this regime, the Floquet bands 



couple and the quasi-energies obtained in Eq|20] are no 
longer valid. Then, we numerically diagonalize Eq J 181 for 
n.m = 15 sidebands, and compare with the calculation 
for the finite system, which requires a larger number of 
Floquet bands in order to reach convergence (n, m = 85). 
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Figure 7: Quasi-energies vs Eo/u) for uj ~ r, r' . We compare the model with periodic boundary conditions (left) and the finite 
one (right). Both models only differ in the finite size effects. The periodic boundary model has considered m,n = 15 and up 
to fifteen order coupling, while the finite size model has considered m,n — 85. Parameters uj — 1.2, A = 0.5, and feo = (all 
parameters in units of r = 1). 



The full Floquet operator in time domain is given by 

(H(k,t) - id t ) \u a>k ) = e a ,k\u a . k )- (21) 

For the calculation of T-L (k,t), the coupling between dif- 
ferent Floquet bands depends on the amplitude of the 
vector potential. If we assume Aq < 1, we can neglect 
the contributions from J p >2 (Aq), leading to a Hamilto- 
nian up to next nearest neighbors: 

Pf (k,t) = A Jo (x) + 2i\ J\ (x) sin (ut) 

+2AJ 2 (x) cos (2ut) + e lkao J Q (y) 

-2e ikao [Ui (y) sin (ut) - J 2 (y) cos (2ut)) , 



which, in the limit Aq — > becomes the energy dispersion 
of the undriven system. 

Eql2T1 is not exactly solvable for the time dependent 
Hamiltonian H (fe, t) NN (Ea l^j) . However, for frequency 
values uj <C t, t' we can neglect the time derivative, such 
that H(k,t) ~ H(k,t) NN . Then, Eqj21] describes an 
eigenvalue equation, and we can classify the mappings 
(k,t) -> U (k,t), from the FFBZ to the set of Floquet 
operators: 

H[k,t)=Tg(k,t)-a. (23) 

The calculation of the first Chern number using Eql23l 
gives ci = 0, because the Hamiltonian belongs to the BDI 
class. However, we can calculate the winding number 
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V\ along the two inequivalent directions of the 2-torus k 
and t as: vx{rf) = 5 $ -^-<f) (//, 77) dfx, where fx, v = k,t, 

(f>{k,t) = arctan^^-^, and g x ^ y are the components of 

the vector g (fc, t) in Eq l23l As we discussed in the main 
text, in this case v\ (i) = 7r for all t, and v\ (k) = for 
all k. 

In the present example of a dimer chain, the existence 



of topologically protected states at the boundary of the 
chain is related with the loops along k, i.e., with the 
winding number v\ (t) . Hence, V\ (t) — n for all t means 
the existence of boundary states for all t, independently 
of the value of Aq , as far as we are in the low frequency 
regime (see FigJS]). 
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Figure 8: (Left) Quasi-energy spectrum within {—u)/2,u/2) for a finite dimers chain in the low frequency regime (uu <C t,t'). 
Right plot shows a zoom of the zero energy mode. The plots show a large number of crossings between the quasi-energies. 
However, the existence of zero energy modes is not affected by these crossings, and they exists for all values of Eq/uj. We have 
considered a finite chain of 20 dimers, n,m— 105, ui = 3.10 -2 , and A = 0.3 in units of r. 



S.3.2: Bands Inversions 

Bands inversions in periodically driven lattices are pro- 
duced due to the coupling between different Floquet 
bands as we decrease the frequency. They change the 
topological properties of the system by opening and clos- 
ing the gap. Fig[9]shows the band structure of the dimers 
chain for different values of uj. The exact crossings be- 
tween conduction and valence band in this model are 
critical points for a topological phase transition, as the 
appearance/disappearance of boundary states demon- 
strate. While the critical points for to = 2 (Fig[9] left) 



at Eq/uj ~ 1.5, and 4.5 are almost fixed as we decrease 
the frequency, the critical point at Eq/lu ~ 2.8 shifts 
to the left as the frequency is lowered (Fig(5] center and 
right). It reaches the origin at uj\ = 1.3 and vanishes 
modifying the phase diagram. Importantly, the shift of 
the critical points depends on the band structure of the 
system under consideration (the critical points can disap- 
pear at the origin, or new ones can show up). However, 
the gap always closes for w„, independently of the critical 
points dynamics (being u n = de/nu, n G Z + , and 8e the 
maximum thickness of the Floquet band). 
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Figure 9: Quasi-energies vs Ao = Eo/ui for different frequencies considering periodic boundary conditions. In dotted green color 
we plot the boundary states, which has been obtained using a finite tight binding model. Left figure shows the quasi-energies 
for uj = 2, where the bands are not inverted. The center figure shows the quasi-energies for the critical value of uj\ — 1.3 at 
which the Floquet bands that belong to close the gap at Ao = 0. The right figure considers a lower frequency (uj = 1) where 
the gap is reopened. Importantly, the comparison between left and right figures shows that the order in which the boundary 
states appear as a function of Ao is inverted. 
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